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ABSTRACT Small RNAs play a critical role in host-pathogen interaction. Indeed, small RNA-mediated silencing or RNA interfer- 
ence (RNAi) is one of the earliest forms of antiviral immunity. Although it represents the main defense system against viruses in 
many organisms, the antiviral role of RNAi has not been clearly proven in higher vertebrates. However, it is well established that 
their response to viral infection relies on the recognition of viral RNAs by host pattern recognition receptors (PRRs) to trigger 
activation of the interferon pathway. In the present work, we report the existence of a novel small noncoding RNA population 
produced in mammalian cells upon RNA virus infection. Using Sindbis virus (SINV) as a prototypic arbovirus model, we pro- 
filed the small RNA population of infected cells in both human and African green monkey cell lines. Here, we provide evidence 
for the presence of discrete small RNAs of viral origin that are not associated with the RNA-induced silencing complex (RISC), 
that are highly expressed and detected by Northern blot analysis, and that accumulate as 21- to 28-nucleotide (nt) species during 
infection. We report that the cellular antiviral endoribonuclease RNase L cleaves the viral genome, producing in turn the small 
RNAs. Surprisingly, we uncovered the presence of a modification on the 3 '-end nucleotide of SINV-derived viral small RNAs 
(SvsRNAs) that might be at the origin of their stability. Altogether, our findings show that stable modified small viral RNAs 
could represent a novel way to modulate host-virus interaction upon SINV infection. 

IMPORTANCE In a continuous arms race, viruses have to deal with host antiviral responses in order to successfully establish an 
infection. In mammalian cells, the host defense mechanism relies on the recognition of viral RNAs, resulting in the activation of 
type I interferons (IFNs). In turn, the expression of many interferon-stimulated genes (ISGs) is induced to inhibit viral replica- 
tion. Here we report that the cytoplasmic, interferon-induced, cellular endoribonuclease RNase L is involved in the accumula- 
tion of a novel small RNA population of viral origin. These small RNAs are produced upon SINV infection of mammalian cells 
and are stabilized by a 3 '-end modification. Altogether, our findings indicate that in our system RNA silencing is not active 
against Sindbis virus (SINV) and might open the way to a better understanding of the antiviral response mediated by a novel 
class of small RNAs. 
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In order to successfully establish an infection, viruses engage in a 
continuous arms race with their hosts, the outcome of which 
might ultimately result in either viral persistence, clearance of the 
viral infection, or death of the infected cells. The arsenal of anti- 
viral defenses and viral counter defense mechanisms is vast and 
very diverse among the different phyla of life. Plants, insects, nem- 
atodes, and fungi mostly (if not exclusively) rely upon an innate 
immune response, whereas mammals have evolved an adaptive 
response in addition to this first line of defense. In the former 
organisms, the antiviral component of the innate immune re- 
sponse is mainly based on a highly conserved process commonly 
referred to as RNA silencing, or RNA interference (RNAi). In this 
phenomenon, long double-stranded RNAs of viral origin are 
cleaved by the RNase Dicer into small interfering RNAs (siRNAs), 
which are then assembled into effector complexes that contain a 



member of the Argonaute (Ago) family; this results in targeted 
degradation of viral messenger or genomic RNAs (1). 

In mammals, another type of innate immunity has evolved to 
control viruses (2). The initial recognition of nonself viral nucleic 
acids by extra- and intracellular sensors triggers the activation of 
type I interferons (IFNs) (3). This leads in turn to the upregulation 
of numerous IFN-stimulated genes (ISGs) (4), which either di- 
rectly trigger a signaling cascade directed against the viral infec- 
tion or mobilize other cells of the immune system. ISGs encode 
proteins involved in apoptosis induction, protein synthesis block- 
ing, or in the regulation of mRNA editing or RNA degradation ( 5 ) . 
One key factor involved in IFN type I-mediated host defense is 
RNase L, a latent cytoplasmic endoribonuclease that is activated 
by 2',5'-oligoadenylates in response to double-stranded RNA 
(dsRNA) sensing and that cleaves both viral and cellular RNAs (6). 



November/December 2013 Volume 4 Issue 6 e00698-13 



Bio' mbio.asm.org 1 



Girardi et al. 



Although extensive deep sequencing analyses have been per- 
formed to identify siRNAs in mammalian cells infected with sev- 
eral RNA viruses, there is so far no good evidence for a possible 
antiviral role of RNAi in mammalian somatic cells (7). Neverthe- 
less, another class of small noncoding RNAs, microRNAs 
(miRNAs), has been shown to play an important role during viral 
infection in mammals. miRNAs are evolutionarily conserved 
small RNAs derived from large primary transcripts, which are 
sequentially processed by the respective nuclear and cytoplasmic 
RNase III enzymes, Drosha and Dicer, to generate mature single- 
stranded -21- or 22-nucleotide (nt) RNAs (8, 9). Similar to 
siRNAs, they are incorporated into an effector Ago-containing 
RNA-induced silencing complex (RISC), whereby they mediate 
posttranscriptional regulation of target mRNAs via partially com- 
plementary sites (lO).Ononehand, miRNAs of cellular origin can 
directly or indirectly regulate viral infections (11, 12). On the 
other hand, some viruses have evolved the capacity to encode their 
own miRNAs, which represent an ideal tool to stealthily modulate 
the cellular environment. To date, viral miRNAs have been almost 
exclusively identified in the genomes of DNA viruses, mostly her- 
pesviruses (13-17), with the notable exception of bovine leukemia 
virus, a retrovirus with an RNA genome (18). It has been com- 
monly assumed that cytoplasmic RNA viruses cannot encode 
miRNAs, not only because they do not have access to the nuclear 
biogenesis machinery but also because their genomic integrity 
would be destabilized by miRNA processing (19). Nevertheless, it 
has recently been demonstrated that the insertion of a cellular 
miRNA precursor into the 3 '-end nontranslated region (NTR) of 
Sindbis virus (SINV) genomic RNA leads to its cytoplasmic pro- 
cessing without affecting the viral replication (20, 21). 

The arthropod-borne SINV is a small, enveloped, positive, 
single-stranded RNA virus and is the prototype for the alphavirus 
genus. Alphaviruses represent a group of widely distributed hu- 
man and animal pathogens, which pose a serious public health 
threat (22). Some of them induce febrile and arthritogenic dis- 
eases, while others can cause highly debilitating diseases, such as 
encephalitis. The SINV genomic RNA is capped and polyadenyl- 
ated and is infectious as naked RNA. Upon entry into the cyto- 
plasm by endocytosis, the host translational machinery recognizes 
the genomic RNA, and four nonstructural proteins are produced 
(23). Their expression is sufficient for the establishment of repli- 
cation complexes and the synthesis of a complementary antigeno- 
mic negative RNA strand. The antigenome is necessary for the 
replication of the viral genome and the production of a sub- 
genomic RNA that ensures the translation of the structural pro- 
teins during the late phase of infection (24, 25). 

To date, SINV has not been shown to give rise to any small 
RNA species in mammalian cells. As such, we decided to investi- 
gate the small RNA profile of SINV-infected human and African 
green monkey cell lines. Using small RNA cloning and deep se- 
quencing techniques, we provide evidence that the SINV genome 
is a source of small RNAs in infected mammalian cells. Some of 
these small RNAs accumulate to levels detected by Northern blot 
analysis. We set out to identify the factors involved in their bio- 
genesis and show that they do not depend on the RNA silencing 
machinery, but rather, that they are downstream products of the 
cytoplasmic IFN-induced endoribonuclease RNase L. Finally, we 
also present experimental evidence that these viral small RNAs are 
modified at their 3' extremity. Altogether, our results indicate that 



the SINV-derived small RNAs could represent key markers of the 
host defense mechanism. 

RESULTS 

SINV-derived small RNAs accumulate in infected mammalian 
cells. In order to investigate the presence and nature of virus- 
derived small RNAs in mammalian cells, we used a small RNA 
cloning and deep sequencing approach. We generated and ana- 
lyzed small RNA libraries from both HEK293 and Vero cells in- 
fected with SINV or not infected with SINV. 

Total RNA was isolated from uninfected cells or cells infected 
with SINV at a multiplicity of infection (MOI) of 0.01 plaque 
forming units (PFU) per cell 16 h postinfection (hpi). Using these 
conditions, we could detect cells with a high viral load without 
observing any apparent cell death as assessed by propidium iodide 
staining (see Fig. SI in the supplemental material). We also con- 
firmed that until 24 hpi, cells were still actively proliferating (data 
not shown). High-throughput sequencing yielded 201,793 and 
1,193,893 reads that mapped to the SINV genome in infected 
HEK293 and Vero cells, respectively (Fig. 1 A); these two numbers 
of reads represented a percentage of 1.66 and 7.99%, respectively, 
of total mapped small RNAs in HEK293 and Vero cells. The 
greater fraction of viral reads in Vero cells is consistent with the 
fact that these cells are known to be more permissive to viral in- 
fections due to a defect in interferon production (26). The distri- 
bution of the reads on the viral genomic (positive) and antigeno- 
mic (negative) strands is indicated in Fig. 1A. Interestingly, the 
vast majority of the reads (98 to 99%) originate from the genomic 
strand of the virus. Nonetheless, we identified limited sequencing 
reads that were derived from the negative strand, albeit at ex- 
tremely low levels in both libraries (Fig. 1A). Their genomic dis- 
tribution showed a peak at the 3' end of the antigenome (see 
Fig. S2A in the supplemental material). This specific region cor- 
responds to the viral promoter, necessary for the replication of the 
viral genome and is known to fold into a conserved 44-nt stem- 
loop structure (25). The length distribution of the antigenomic 
small RNA reads peaked at a size of 22 nt (Fig. S2B), which corre- 
sponds to the 5' arm of the stem-loop structure (Fig. S2C). Al- 
though we could validate the presence of the 44-nt-long RNA by 
Northern blot analysis, we were unable to detect the accumulation 
of the small 22-nt RNA as a discrete species (Fig. S2D). 

The size distribution of the main population of viral reads 
mapping on the positive strand of the virus was broad with no 
defined peak at a specific length (Fig. IB). However, a closer in- 
spection revealed that the overall size and genomic distribution of 
the viral reads were highly similar in both HEK293 and Vero cells 
(Fig. 1C). Moreover, a greater number of reads were found be- 
tween nucleotides 7000 and 11700 of the viral genome, corre- 
sponding to the abundantly expressed subgenomic RNA. 

We selected 14 candidates among the most abundant peaks 
obtained in both cell types for validation by Northern blot analysis 
(see Fig. S3A in the supplemental material). We could detect sev- 
eral SINV-derived viral small RNAs (SvsRNAs) of 20 to 30 nt in 
length that were derived from all along the viral genome 
(Fig. S3B). Among them, SvsRNA-1, -2, and -3 were detected in 
both HEK293 and Vero cells (Fig. S3C). The presence of abundant 
larger bands accumulating between 40 and 100 nt was also evi- 
denced. Of note, the expression of both larger and smaller viral 
RNAs increased proportionally to the accumulation of the viral 
genome during a time course of infection (Fig. ID and E). The 
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FIG 1 Identification of SINV-derived small RNAs. (A) Virus-specific small RNA reads were recovered from each library (HEK293 cells infected with SINV 
[HEK293-SINV] and Vero-SINV) and were mapped to both positive ( + ) and negative ( — ) strands of the viral genome. (B) Size distribution of SINV-derived 
small RNA populations. The proportion of viral sRNAs that were derived from the genomic sense strand in each size class is shown as a percentage of the total 
viral small RNA reads for this strand. (C) Distribution of viral reads mapping to the SINV genomic strand. The abundance of small RNAs was calculated and 
plotted as the sum of normalized reads (per 10 5 viral reads) in each single-nucleotide sliding window along the viral genome. The small black arrows indicate the 
peaks chosen for Northern blot validation (see Fig. SI in the supplemental material). A schematic diagram represents the organization of SINV genome. The two 
open reading frames (ORF), which encode the nonstructural (ns) and structural proteins, are shown. The nontranslated region (NTR) at the 3' end of the virus 
is shown as a small black bar. (D) Northern blot validation of SINV candidates. Total RNA was isolated at the indicated time points from HEK293 cells infected 
with SINV at an MOI of 0.01. Noninfected cells ( — ) were used as a negative control. U6 snRNA was used as a loading control. (E) Quantification of the SINV 
genome during infection by reverse transcription-quantitative PCR (RT-qPCR). The relative quantification is normalized to glyceraldehyde-3-phosphate 
dehydrogenase (GAPDH) levels and is represented on a logarithmic scale as a mean of three replicates. NI, not infected. 



identification of SINV-derived small RNAs prompted us to ana- 
lyze more closely their origin and biochemical properties. 

The production of SvsRNAs does not rely on the miRNA bio- 
genesis machinery. The identification of small RNAs from RNA 
viruses could be seen either as a result of a productive mechanism 
for the pathogen (as for virus-encoded miRNAs) or as the conse- 
quence of a host response to control the infection by degrading 
viral RNAs. SINV expresses four enzymes that play key functions 
during the viral replicative cycle. However, none of these enzymes 
is known to possess an endoribonuclease activity (22) that could 
enable the processing of small RNAs from a viral RNA precursor. 
Therefore, cellular enzymes have to be involved in the biogenesis 
of SINV sRNAs. In order to identify them, we decided to knock 
down cytoplasmic enzymes with a known RNase activity. We first 
focused on proteins important for the production of small RNAs, 
such as miRNAs. Hence, we transfected siRNAs against Drosha, 
Dicer, or Ago2 followed by an infection with SINV and an analysis 
of viral small RNA accumulation. Although Drosha is known to be 
mostly nuclear, we also took it into consideration because it was 
recently shown to relocalize to the cytoplasm upon SINV infection 
(21). Although we were able to downregulate very efficiently 
Drosha, Dicer, and Ago2 (see Fig. S4A to C in the supplemental 
material), the knockdown of these factors neither affected the viral 



load nor the accumulation of SvsRNA-1, -2 and -3 (Fig. S4D to G). 
We also examined whether we could detect an siRNA signature 
among all viral reads by looking for small RNAs that could per- 
fectly base pair on 19 or 20 nt and present a 2-nt overhang in 3'. A 
very limited number of such duplexes could be found, but they 
were not more represented than larger duplexes presenting with 3 ' 
or 5' overhangs (data not shown). It has been reported that small 
RNAs, which are not generated by Dicer or Drosha cleavage, could 
nevertheless be assembled into Argonaute proteins (e.g., tRNA 
fragments [27]). We therefore immunoprecipitated Argonaute 
proteins followed by Northern blot analysis to see whether 
SvsRNAs could be detected within the RISC. However, although 
we could detect the accumulation of the cellular miRNA miR-16, 
using both an antibody specific for Ago2 (Fig. S4H, left panel) and 
an antibody able to recognize all four Ago proteins (Fig. S4H, right 
panel) , we were unable to detect SvsRNAs within immunoprecipi- 
tated RISC. 

RNase L is involved in the production of SvsRNAs. We then 
hypothesized that an endoribonuclease involved in the antiviral 
response could be implicated in SvsRNA production. RNase L is a 
ubiquitous, cytoplasmic, interferon-induced endonuclease that 
plays a key role in the cellular response to viral infections. It cleaves 
single-stranded RNA regions after the dinucleotides UA/UU, 
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FIG 2 Effects of RNase L and Xrnl knockdown on SvsRNA accumulation. (A 
and D) RNase L and Xrnl knockdown in SINV-infected HEK293 cells was 
verified by Western blotting. Tubulin (Tub) was used as a normalizer. sictrl, 
control siRNA; siRNaseL, RNase L-specific siRNA; siXrnl, Xrnl-specific 
siRNA. (B and E) SINV genomic RNA was quantified by RT-qPCR, before and 
after RNase L or Xrnl KD. The relative quantification is normalized to 
GAPDH levels and is represented on a logarithmic scale as a mean of three 
replicates. (C and F) Northern blot analysis on total RNA from infected 
HEK293 cells before and after RNase L or Xrnl KD. U6 was used as a normal - 



thereby preventing virus accumulation (28). As such, we knocked 
down RNase L using siRNAs, before infecting the cells with SINV. 
As expected, RNase L knockdown (KD) resulted in the upregula- 
tion of the SINV genome by approximately 2.5-fold (Fig. 2A and 
B). Interestingly, it also resulted in a dramatic reduction of the 
viral small RNAs SvsRNA- 1, -2, and -3 and their precursors 
(Fig. 2C). 

In order to extend this observation and assess the role of RNase 
L in SvsRNA generation, we also produced small RNA libraries 
from HEK293 cells treated with a control or an RNase L-specific 
siRNA prior to SINV infection. In accordance with the striking 
effect seen on the three specific SvsRNAs, reduction in RNase L 
protein level decreased the entire population of clonable viral se- 
quences from 2.75 to 1.19% (see Fig. S5A in the supplemental 
material). However, RNase L knockdown affected neither the size 
distribution of viral reads nor their genomic distribution (Fig. S5B 
and C). 

Given the cleavage properties of RNase L (see next section), it 



was unlikely that the identified small RNAs were direct products 
of this enzyme. We therefore looked for other factors implicated 
in RNA degradation that could affect SvsRNA accumulation. We 
used RNAi to downregulate the expression of the major cytoplas- 
mic5'-3' exoribonucleaseXrnl (Fig. 2D). As expected, thecapped 
viral genome was not affected by downregulation of Xrnl 
(Fig. 2E). However, the levels of SvsRNA-1, -2, and -3 and their 
precursors were stabilized (Fig. 2F) . We then looked at the effect of 
combining knockdown of these factors (see Fig. S6 in the supple- 
mental material). We found that the double knockdown of both 
RNase L and Dicer did not affect the accumulation of SvsRNA- 1 
compared to the single RNase L knockdown, which confirms that 
Dicer does not play a role (Fig. S6A). On the contrary, the accu- 
mulation of SvsRNA-1 upon double knockdown of both RNase L 
and Xrnl increased compared to the single RNase L knockdown 
(Fig. S6A). This finding suggests that both RNase L and Xrnl can 
act upon viral RNAs and have opposite activities in terms of gen- 
eration and stabilization of the SvsRNAs. Xrnl might act up- 
stream of RNase L in degrading, for example, uncapped viral 
RNAs that would otherwise be a substrate for RNase L; alterna- 
tively, it could directly act upon the RNase L products. The former 
hypothesis is more likely though because Xrn 1 has a preference for 
5' phosphorylated RNAs, whereas RNase L leaves a 5' OH after 
cleavage. We also performed a high-molecular-weight Northern 
blot analysis to look at the effects of the various knockdowns on 
the accumulation of the genomic and subgenomic viral RNAs. 
The results confirm the quantitative PCR data, in that only the 
RNase L KD has an impact on the accumulation of both genomic 
and subgenomic RNAs (Fig. S6B). 

We also tested several cytoplasmic 3 '-5' exoribonucleases, in- 
cluding DIS3L, DIS3L2, and Rrp40 (also known as EXOSC3, a 
core sub unit of the exosome); knockdown of these factors had no 
effect on SvsRNA accumulation (data not shown). This indicated 
that at least the main cytoplasmic 3 '-5' exoribonucleases do not 
seem to be implicated in the production of the viral small RNAs. 

The 3' extremity of SvsRNA-1 is modified. RNase L process- 
ing produces a typical signature: the enzyme leaves a 5' hydroxyl 
(OH) and a 3' phosphate (P) after cleavage (6). However, these 
features are not compatible with our small RNA cloning protocol 
(29), which strictly relies on the presence of a 5' P and a 3' OH on 
the small RNA extremities. This prompted us to characterize the 
extremities of SvsRNAs to better understand their biogenesis. We 
initially decided to restrict our analysis to SvsRNA- 1 because it 
accumulated as a major form of 23 nt. The assay used for the 
analysis of RNA 5' ends involved digestion with the Terminator 
exonuclease, a processive 5'-3' exonuclease that specifically tar- 
gets and degrades RNAs with a 5' monophosphate, but not RNAs 
with 5' triphosphates, 5' OH, or a 5' CAP (30). 

For a control, the Terminator reaction was carried out on both 
a 5' OH and a 5' P oligoribonucleotide with the same sequence as 
SvsRNA- 1. Terminator treatment did not result in the degrada- 
tion of the 5' OH synthetic oligoribonucleotide, whereas the 5' P 
RNA oligonucleotide was almost completely degraded. Accord- 
ingly, SvsRNA- 1 is prone to degradation by the Terminator exo- 
nuclease, which is consistent with the presence of a monophos- 
phate at the 5' end, while the accumulation of the precursor bands 
was less affected by the treatment (Fig. 3A). 

Next, we performed sodium periodate (NaI0 4 ) treatment 
followed by j3-elimination to specifically examine the 3' end of 
SvsRNA- 1. The presence of hydroxyl groups in positions 2' and 
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FIG 3 Characterization of 5' and 3' extremities of SvsRNA-1. (A) Synthetic 
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3' on the last ribose is needed for the periodate reaction. 
j8-elimination of the oxidized RNA results in an RNA product one 
or two nucleotides shorter. On the contrary, the presence of mod- 
ifications on the ribose of the 3' terminal residue renders the RNA 
insensitive to periodate oxidation and impairs the chain scission 
(31). For a control for the reaction, we measured its effect on the 
unmodified cellular miRNA miR- 1 6. As expected, the electropho- 
retic mobility of the miRNA strongly changed after treatment 
(Fig. 3B). In contrast, SvsRNA- 1 mobility was found to be insen- 
sitive to periodate oxidation and j3-elimination (Fig. 3B). This 
indicated the presence of a modification on the terminal ribose, 
most likely in the 2' position, which is perfectly compatible with 
our cloning protocol. Interestingly, it also seems that the interme- 
diate RNA products are insensitive to periodate oxidation and are 
consequently modified (see Fig. S7 in the supplemental material). 
Moreover, we were also able to assess the existence of a modified 3 ' 
end for other SvsRNAs by Northern blotting (Fig. S8). 

In order to gain more insight into the nature of this modifica- 
tion, we treated our samples with calf intestinal phosphatase (CIP) 
(Fig. 3C, lanes 2, 4, 6, 8, 10, 12, 14, and 16) prior to jS- elimination 
(Fig. 3C, lanes 4, 8, 12, and 16). The double treatment strongly 
affected the mobility of a control oligoribonucleotide with a 3' P 
(Fig. 3C, lane 4), while the mobility of the 3' OH oligonucleotide 
was affected independently of the CIP treatment (Fig. 3C, lanes 7 
and 8). Of note, the combined reactions changed neither the mo- 
bility of a control oligoribonucleotide bearing a CIP-insensitive 



modification (methyl group in position 2' of the last ribose) 
(Fig. 3C,lane 12) nor the mobility of SvsRNA- 1 (Fig. 3C,lane 16). 

Taken together, these results show that SvsRNA- 1 is modified 
at its 3' extremity but that the modification is not due to a mono- 
phosphate group caused by a direct RNase L cut. 

Viral genomic RNA is modified on A 8333 . The presence of 
posttranscriptional modifications, such as 2'-0-methylation, is a 
characteristic of several cellular RNAs in mammals, such as rRNAs 
(32), snRNAs (33), and Piwi-interacting RNAs (piRNAs) (34, 35). 
Recently, it has also been reported for several viruses that their 
genomic RNA can be 2'-0-methylated at specific positions (36, 
37). We hypothesized that the modification on the small RNA 
could be derived from a modification at the level of the genomic 
RNA, which is reasonable given the fact that SvsRNA- 1 precursor 
RNAs also appear to be modified (see Fig. S7 in the supplemental 
material). We thus looked for modifications in the genomic re- 
gion spanning SvsRNA- 1 using the RTL-P (reverse transcription 
at low deoxyribonucleoside triphosphate concentrations followed 
by polymerase chain reaction) technique (38): the tendency of the 
reverse transcriptase to pause immediately before a 2'-0- 
methylated nucleotide is exploited to establish the efficiency of 
retrotranscription by semiquantitative PCR, using differential de- 
oxynucleoside triphosphate (dNTP) concentrations. 

As shown in Fig. 4A, we designed a reverse transcription 
primer specific for the viral genome, as well as primers for semi- 
quantitative PCR in the region surrounding SvsRNA- 1 (shown in 
red): one reverse primer (R) and three forward primers, upstream 
(F v ) or downstream of the putative modification (F D1 and F D2 ). 
After retrotranscription with low and high dNTP concentrations, 
we performed multiplex PCR using primers Fu/Fd (in the 
SvsRNA- 1 region) or primers F D1 /F D2 (in the control region) and 
then compared the signal intensity ratios of PCR products ampli- 
fied with different numbers of PCR cycles (Fig. 4B). Our data 
indicated that only the Fjj/Fj-,) ratio changed depending on the RT 
conditions (Fig. 4B, top panel), suggesting that at a low dNTP 
concentration, the reverse transcriptase paused in the SvsRNA- 1 
region, most probably due to the presence of a modification on the 
viral genome in the region of interest. 

To further strengthen this observation and to identify the pre- 
cise position of the nucleotide modification, we performed primer 
extension assays on RNA isolated from either SINV-infected 
HEK293 or Vero cells (Fig. 4C). 2'-0-methylation can cause re- 
verse transcriptase to pause one nucleotide before and/or at the 
O-methylated nucleotide. Indeed, we identified two bands (RT 
stops) that correspond to nucleotides G 8334 and A 8333 , the latter 
corresponding precisely to the 3' end of the SvsRNA- 1 sequence 
(Fig. 4C). 

These data imply that the 3 '-end modification of SvsRNA- 1 
was likely already present on the viral genome and might be a 
2 ' - O-methylation. 

3 '-end modification is a hallmark of SvsRNAs. Taken to- 
gether, our observations indicate that the SINV genomic RNA is 
modified, and are therefore suggestive of a potential modification 
of SvsRNAs. We then decided to investigate to what extent reads 
mapping to the SINV genome were modified at the last nucleotide 
position. Small RNA libraries were prepared from noninfected or 
SINV-infected HEK293 cells using total RNA treated (or not) with 
NaI0 4 . This treatment specifically impairs the 3' adapter ligation 
to RNA species with unmodified 3' ends (2' and 3' OH), while 
RNAs bearing a modification at one of these positions are resistant 
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to oxidation, and can therefore be ligated. Overnight ligation of 
the 3' adapter at 16°C was performed in order to enhance the 
cloning of 2'-0 modified small RNAs (39). 

After deep sequencing of these libraries, we mapped the total 
number of reads against both the human and SINV genome (Ta- 
ble 1). As expected, in both noninfected and infected cells, the 
NaI0 4 treatment drastically reduced the number of clonable 
miRNA sequences (see also Fig. S9A in the supplemental mate- 
rial). In contrast, modified small RNAs (such as tRNA fragments) 
were preferentially cloned over miRNAs upon NaI0 4 treatment 
(Fig. S9B). Strikingly, almost 75% of the viral reads were kept after 
oxidation (Table 1 ) . The oxidation treatment did not dramatically 
change the genomic distribution of the viral reads compared to the 
nontreated sample (Fig. 5A). However, deep sequencing revealed 
some new peaks (Fig. 5A, arrows) that we could also validate by 
j8-elimination and Northern blot analysis (see SvsRNA- 15 and - 16 
in Fig. S8). 

The size distribution of the viral reads remained unchanged 
before and after NaI0 4 treatment, indicating that SINV-derived 
modified RNAs belonged to each length (Fig. 5B). Finally, we an- 
alyzed the nucleotide composition of the 3' end of SINV small 
RNAs and observed that we had enrichment in adenosine (A) at 
the 3' end after oxidation, which was not clearly discernible in the 
nontreated sample (Fig. 5C). This is in line with the specific 2'-0- 
methylation occurring on internal adenosines recently identified 
on the genome of dengue virus (36). 
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FIG 4 SINV genomic RNA is modified. (A) Schematic diagram of the primer 
design for RTL-P to detect the presence of 2'-0-methylated nucleotides (m) 
on the viral genome. SvsRNA- 1 is indicated in red. (B) Detection of methyl- 
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DISCUSSION 

Mammals have evolved a very sophisticated innate immune re- 
sponse to ward off pathogens that is based on the sensing of 
pathogen-associated molecular patterns (PAMPs) by dedicated 
pattern recognition receptors (PRRs). When it comes to viral in- 
fection, one of the most potent PAMPs is dsRNA. The recognition 
of dsRNA by PRRs triggers a signaling cascade that ultimately 
leads to the activation of interferon-sensitive genes (ISGs) (40). 
RNA silencing is another type of antiviral innate immunity, 
which, as of now, has mostly been demonstrated in plant and 
insect organisms (1). Nonetheless, a role for RNA silencing as a 
host defense system against mammalian viruses could reasonably 
be considered. Indeed, Dicer, the key enzyme for dsRNA cleavage 
is conserved in higher vertebrates, and although its main substrate 
in vivo is pre-miRNA molecules, it is able to process long dsRNA 
into short RNAs in vitro (41). However, the question as to which 
other types of dsRNA molecules mammalian Dicer can cleave in 
vivo has been the matter of intense research, and especially so 
when considering RNAs of viral origin. Here, we wanted to make 
our contribution to the field by studying the prototypical arbovi- 
rus, Sindbis virus. This type of virus can infect both insects and 
mammals, and as such is a useful tool to study the conservation of 
antiviral mechanisms in both phyla. Intriguingly, it has been 
shown for some members of the alphavirus genus that their 



Figure Legend Continued 

ation by RTL-P under both high and low dNTP concentrations and PCR cycles 
on SINV-infected HEK293 cells. The ratio of PCR signal intensity is indicated 
beneath each lane. (C) Primer extension mapping of 2'-0-methylated nucle- 
otides. Lanes U, G, C, and A represent dideoxy sequencing reactions per- 
formed on the in vitro transcript amplified for the same viral region. Asterisks 
indicate stops of the reverse transcriptase. 
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TABLE 1 Deep sequencing results of Nal0 4 -treated and nontreated small RNA libraries" 



Category'' 


HEK293-NI-NT 


HEK293-NI-NalC 


'4 


HEK293-SINV-NT 


HEK293-SINV-Nal0 4 


Read count 


Read % 


Read count 


Read % 


Read count 


Read % 


Read count 


Read % 


Raw data 


21,631,645 




34,131,455 




20,753,238 




28,263,345 




Filtered out 


92,556 




3,589,935 




146,674 




5,425,669 




Unmapped 


840,747 




7,135,376 




1,276,025 




5,143,518 




Mapped 


20,698,342 


100.00 


23,406,144 


100.00 


19,330,539 


100.00 


17,694,158 


100.00 


hsa miRNAs 


18,298,024 


88.40 


7,325,147 


31.30 


14,850,144 


76.82 


4,138,513 


23.39 


hsa rRNAs 


66,614 


0.32 


1,509,137 


6.45 


337,226 


1.74 


850,630 


4.81 


hsa tRNAs 


878,642 


4.24 


6,260,105 


26.75 


1,680,819 


8.70 


5,915,996 


33.43 


hsa snRNAs 


8,956 


0.04 


13,719 


0.06 


8,913 


0.05 


6,996 


0.04 


hsa snoRNAs 


91,453 


0.44 


93,919 


0.40 


111,048 


0.57 


70,898 


0.40 


hsa miscRNAs 


59,911 


0.29 


36,831 


0.16 


283,636 


1.47 


98,472 


0.56 


Only viral RNA 


491 


0.00 


37 


0.00 


400,252 


2.07 


262,916 


1.49 


Unannotated RNA 


1,294,251 


6.25 


8,167,249 


34.89 


1,658,501 


8.58 


6,349,737 


35.89 



" HEK293 cells were infected with SINV or not infected (NI) and treated with sodium periodate (Nal0 4 ) or not treated with sodium periodate (NT). 
b Abbreviations: hsa, homo sapiens; miscRNAs, miscellaneous RNAs. 



genomic RNA can be processed into both siRNAs and piRNAs in 
insect somatic cells (42, 43). 

In this study, we generated small RNA libraries from two dif- 
ferent mammalian cell lines infected with SINV. Although the 
percentage of sequences mapping to the viral genome was signif- 
icant (-1.7 to ~8% in HEK293 and Vero cells, respectively), we 
were unable to detect either siRNA or piRNA signatures. On the 
contrary, the viral reads displayed a highly biased strand distribu- 
tion. Although antigenomic (negative) strand-derived sequences 
represented a very small percentage of the total viral reads, they 
showed properties distinct from the genomic (positive) strand 
reads: their size distribution peaked at 22 nt, and their localization 
on the viral genome is restricted to the promoters required for the 
synthesis of both the genomic and subgenomic RNA. Although 
the cloned 22-nt sequences were not visible by Northern blot anal- 
ysis, we could detect the accumulation of an ~44-nt product that 
in our knowledge has never been observed before. We hypothesize 
that these RNA species could be critical for viral replication, sim- 
ilarly to what has been previously reported for influenza virus 
(44). 

Consistent with the relative abundance of the positive strand 
compared to the negative strand, the majority of reads were de- 
rived from the plus-strand RNA and more specifically from the 
region corresponding to the subgenomic RNA sequence. We 
could detect several hot spots of small RNA accumulation, indi- 
cating that discrete SINV-derived viral small RNA species (SvsR- 
NAs) could be produced in infected mammalian cells. These Svs- 
RNAs do not display any peak in size distribution and do not map 
to any alphavirus conserved regions. Nonetheless, we noticed that 
they originate from the same viral genomic positions in both in- 
fected African green monkey and human cells. In order to confirm 
that this was not solely due to a bias in the small RNA cloning 
protocol, we performed Northern blot analyses of several small 
RNA candidates. For some of the candidates, we were able to de- 
tect discrete bands, which increased in intensity throughout the 
time course of infection. 

We therefore focused on the characterization of the positive- 
strand-derived small RNAs. One hypothesis regarding their iden- 
tity could be that some of these SvsRNAs are in fact virus-encoded 
miRNAs. Using bioinformatics tools such as miRDeep2 (45, 46) 



and miRanalyzer2 (47, 48), we excluded the possibility that puta- 
tive miRNA precursor structures were present in the viral genome. 
Nonetheless, we tested whether the miRNA machinery was re- 
quired for the production of viral small RNAs by knocking down 
Drosha, Dicer, and Ago2. The downregulation of these factors 
affected neither SINV replication nor accumulation of the three 
most abundant identified SvsRNAs, indicating that their produc- 
tion must involve another RNase. We discovered that the RNase L 
endonuclease is implicated in SvsRNA genesis. Indeed, its down- 
regulation causes a dramatic reduction in the accumulation of 
Sindbis viral small RNAs detected by Northern blot analysis; this is 
accompanied by an upregulation of the viral genomic RNA. In 
order to extend this observation, we also generated small RNA 
libraries from HEK293 cells treated with siRNAs against RNase L 
and infected with SINV. We found that, under these conditions, 
the global population of SINV sRNAs was downregulated. RNase 
L is an endoribonuclease that is present in the cell in the form of an 
inactive monomer. Upon type I interferon induction, and more 
specifically 2'-5' oligoadenylate synthetase activation (6), it can 
dimerize and then act upon its viral RNA targets. RNase L is 
known to play an antiviral role against viruses such as West Nile 
virus (49), hepatitis C virus (50), and SINV (5 1 ). We also looked at 
the involvement of exoribonucleases in the production of 
SvsRNAs and found that, indeed, the accumulation of these small 
RNAs is counteracted by the cellular 5 ' -3 ' exonuclease Xrn 1 . This 
suggests that SvsRNAs represent stable viral products of cellular 
decay enzymes. We therefore hypothesize that RNase L is the en- 
zyme primarily involved in the generation of small RNAs and that 
the further processing operated by other uncharacterized enzymes 
could give rise to the final 20- to 30-nt-long viral RNAs. It also 
remains possible that some of the fragments we detected are direct 
products of RNase L but are further modified by kinases and de- 
phosphorylases, which would allow their cloning. 

The fact that RNA fragments generated by RNase L are de- 
tected implies that some of them must be stabilized. Indeed, we 
discovered that the viral small RNA candidates appear to be mod- 
ified on the 3' proximal nucleotide. We could extend this pecu- 
liarity to the global small RNA population by cloning and se- 
quencing small RNAs after an oxidation treatment designed to 
prevent ligation of nonmodified RNAs. Our data indicate that this 
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FIG 5 Deep sequencing upon periodate oxidation in SINV-infected HEK293 
cells. (A) Distribution of viral reads mapping to the SINV genome. The abun- 
dance of small RNAs was calculated and plotted as the sum of normalized reads 
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RNA. The small black arrows indicate the peaks chosen for Northern blot 
validation. (B) Size distribution of SINV-derived small RNA populations. The 
percentage of viral sRNAs deriving from the genomic sense strand in each size 
class is shown as a percentage of the total viral small RNA reads for this strand. 
Data from SINV-infected HEK293 cell samples treated with NaI0 4 or not 
treated with NaI0 4 (NT) are shown. (C) Relative nucleotide frequency per 
position. The viral reads with 0 mismatch to the SINV positive strand were 
considered for this analysis. Only the last five nucleotides from the 3' end of 
each sequence are displayed. 



modification is most likely a 2 ' - O-methylation and seems to occur 
at the level of the viral genomic RNA, preferentially on adenosines. 
Thanks to a double approach of RTL-P and primer extension, we 
were able to map such a modification on position A 8333 of the 
SINV viral genome. Moreover, our sequencing data indicate that 
there are most likely many others given the enrichment in A resi- 
dues at the 3' end of small RNAs cloned after oxidative treatment. 

The main protective mechanism against 3 '-5' degradation is a 
2' -O-methylation on the 3' proximal ribose, which prevents deg- 
radation by exoribonucleases and thus helps to stabilize small 
RNAs such as miRNAs and siRNAs in plants, piRNAs in animals, 
and siRNAs in Drosophila melanogaster (52). This suggests that 
RNA modification could have evolved as a response to protect the 
anti-pathogen small RNAs during infections. Posttranscriptional 
RNA modifications represent an important but not yet fully un- 
derstood system to regulate transcripts. Cellular RNAs contain 
many distinct posttranscriptional modifications at thousands of 
sites. In the same way, RNA viruses take advantage of nucleoside 
modifications on the base or on the ribose of their genomic RNA 
(36, 53) as a host mimicry strategy to abrogate innate immune 
signaling (54). 

The enzyme responsible for the methylation of SINV RNA re- 
mains unknown. However, it has been reported that the alphavi- 
rus nsP2 protein contains a putative methyltransferase (MT)-like 
domain, similar to dengue virus NS5 MTase (55). Although the 
activity of this domain has not been verified yet, nsP2-MTase mu- 
tant viruses exhibit a strong phenotype: negative-strand synthesis 
is continuous, and unstable viral replication complexes are pro- 
duced. Surprisingly, RNase L-deficient cells infected with SINV 
display the same phenotype, suggesting that the two proteins 
function in similar events (55). For all these reasons, we propose 
that RNase L-dependent, 3 '-modified small RNAs from SINV 
could be involved as intermediates in this mechanism. 

Altogether, our findings show a link between the host degrada- 
tion pathway and stabilization of modified viral small RNAs. 
Clearly, not all the small RNAs identified by our deep sequencing 
approach correspond to RNase L downstream products. None- 
theless, the fact that at least some of these small RNAs are pro- 
tected and detected might implicate them in signaling events, 
maybe upstream of the infection front, which might prove diffi- 
cult to assess in cell culture. We hypothesized that by removing 
RNase L we might uncover signatures corresponding to other 
types of small RNAs, such as siRNAs or piRNAs, but it did not turn 
out to be the case. This could be due to the fact that other enzymes 
might also be involved in the production of RNA fragments and is 
an indication that it will make the detection of antiviral RNAi 
complicated in differentiated mammalian cells. An alternate hy- 
pothesis could also be that SINV encodes an RNA silencing sup- 
pressor that blocks Dicer activity but has no effect on RNase L. 
Finally, it is also possible that the antiviral RNAi response might be 
detected only in specific subtypes of cells in mammals, which 
would be different from the differentiated somatic cells we used 
for our study. Further investigation will be needed to answer these 
questions and to decipher the roles played by this type of virus- 
derived small RNAs, as well as their potential implication in the 
design of new antiviral therapeutic strategies. 

MATERIALS AND METHODS 

Viral stocks, cell culture, and virus infection. Plasmid carrying wild-type 
(WT) and green fluorescent protein (GFP)-SINV genomic sequence 
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(kindly provided by M. C. Saleh) was linearized with Xhol and used as a 
substrate for in vitro transcription using mMESSAGE mMACHINE 
capped RNA transcription kit (Ambion) following the manufacturer's 
instructions. Sindbis viral stock was prepared in BHK21 hamster kidney 
cells, and titers were measured by plaque assay. HEK293 and Vero cells 
were maintained in Dulbecco's modified Eagle medium (DMEM) 
(Gibco) supplemented with 10% fetal bovine serum (FBS) (Clontech) in a 
humidified atmosphere of 5% C0 2 at 37°C. Cells were infected with SINV 
at an MOI of 0.01. Samples were harvested at 0, 4, 8, 12, 16, 24, and 48 h 
postinfection (hpi). 

RNAi-mediated protein depletion. For siRNA transfection, 0.4 X 10 5 
HEK293 cells were seeded in 6-well plates. Transfections were carried out 
for 72 h using 20 nM ON-TARGET plus Smart Pool siRNAs (Dharma- 
con) for one-step depletion (human RNase L and Xrnl) or 100 nM ON- 
TARGET plus Smart Pool siRNAs (Dharmacon) for two-step depletion 
(human Drosha, Dicer, and Ago2) with Lipofectamine 2000 (Invitrogen) 
following the manufacturer's instructions. The ON-TARGET plus Non- 
targeting Pool siRNAs (Dharmacon) was used as a negative control. Cells 
were infected with SINV and harvested after 16 h for RNA and protein 
analysis. 

Small RNA cloning and sequencing. RNA was extracted from nonin- 
fected and SINV-infected HEK293 and Vero cells 16 hpi. Small RNA 
cloning was conducted with 10 or 20 /xg of total RNA as previously de- 
scribed (14, 29). Small RNA libraries from SINV-infected HEK293 cells 
treated with control or RNase L-specific siRNAs were generated following 
the same protocol with some modifications. Mainly, small RNAs were 
ligated using degenerate 5' and 3' adapters to limit biases in the ligation 
step (56). This procedure was also used to prepare small RNA libraries 
after NaI0 4 treatment of 100 /ng of total RNA, as previously described 
(57), except that the 3' ligation was performed at 16°C overnight, in 12% 
polyethylene glycol 8000 (PEG 8000) to favor the cloning of modified 
small RNAs. Sequencing was performed at the IGBMC Microarray and 
Sequencing platform, Illkirch, France, using different Illumina instru- 
ments (Genome Analyzer IIx, HiSeq 2000, and HiSeq 2500) with a read 
length of 36 or 50 bp. 

A detailed description of the deep sequencing data analysis is available 
in Text S 1 in the supplemental material. 

Microarray data accession number. The data discussed in this publi- 
cation have been deposited in NCBI's Gene Expression Omnibus (58) and 
are accessible through GEO Series accession number GSE48831 (http: 
/ /www.ncbi.nlm.nih.gov/geo/ query/acc.cgi?acc = GSE4883 1 ) . 

SUPPLEMENTAL MATERIAL 

Supplemental material for this article may be found at http://mbio.asm.org 
/lookup/suppl/doi: 10.11 28/mBio.00698- 1 3/-/DCSupplemental. 

Text SI, DOC file, 0.1MB. 

Figure SI, TIF file, 1MB. 
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